AD  -  A 1 32  822  FREE  BOUNDARY  PROBLEM  FOR  UNSTEADY  SLAG  FLOW  IN  THE  1/1 

HEARTH(U)  WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH 
CENTER  H  KAWARADA  AUG  83  MRC-TSR-2551  DAAG29-80-C-004 1 
UNCLASSIFIED  F/G  12/1  NL 


I 


END 

10  83 

OTIC 


MICROCOPY  RLSOIU1ION  1E.S1  CHARI 

NA1IONAI  RUHI4U  0»  stAM'AKl's  JUj.  •  A 


ref  CFSV'QV 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 


August  19Ht 


(Received  July  12,  19H3) 


BUS  FILE  COPY 


Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 

l).  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


UNIVERSITY  OP  WISCONS IN-MADISON 
MATHEMATICS  RESEARCH  CENTER 


FREE  BOUNDARY  PROBLEM  FOR  UNSTEADY  SLAG  FLOW  IN  THE  HEARTH 

Hideo  Kawarada* 

Technical  Summary  Report  #2551 
August  1983 


Drainage  of  the  hearth  is  an  important  component  of  a  successful 
operation  of  a  blast  furnace.  In  order  to  investigate  the  influence  of 
tapping  conditions  due  to  the  shape  of  the  slag  surface,  the  one-phase  flow  of 
slag  during  tapping  was  modeled  as  a  free  boundary  problem.  This  problem  was 
reformulated  by  using  the  method  of  integrated  penalty  and,  then  was  simulated 
by  using  a  finite  difference  method  developed  by  the  author. 

The  objective  of  this  report  is  to  give  a  mathematical  justification  of 
the  penalty  method  formulation,  in  which  the  perturbation  with  respect  to  the 
domain  and  the  asymptotic  properties  of  solutions  of  the  boundary  value 
problem  for  an  elliptic  equation  with  penalty  terms  are  used. 
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SIGNIFICANCE  AND  EXPLANATION 


We  can  observe  many  phenomena  involving  Free  Boundaries  in  various  fields 
of  engineering  and  applied  science,  for  example,  jet  problems,  transient 
multi-fluid  flows,  the  equilibria  of  plasmas,  Stefan  problems,  free  boundary 
problems  in  optimal  shape  designs  and  others.  Hence  there  is  interest  and 
need  to  develop  efficient  and  accurate  numerical  methods  for  the  solution  of 
these  problems.  Some  free  boundary  problems  mentioned  above  have  been 
successfully  solved  by  using  the  penalty  method  developed  by  the  author  and 
his  colleagues.  An  important  feature  of  our  approach  is  that  the  outward 
normal  derivative  of  the  solution  at  the  free  boundary  is  approximated 
efficiently. 

Hie  objective  of  this  report  is  to  give  the  required  mathematical 
justification  for  the  model  problem  reformulated  by  the  method  of  integrated 
penalty.  Specifically,  we  prove  the  convergence  of  the  penalised  free 
boundary  to  the  original  one  as  the  penalising  parameter  e  tends  to  sero. 
Here  the  perturbation  theory  with  respect  to  the  domain  and  the  asymptotic 
properties  of  the  solutions  of  a  boundary  value  problem  for  an  elliptic 
equation  with  penalty  terms  are  used.  As  an  application,  the  unsteady  slag 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
suanary  lies  with  HRC,  and  not  with  the  author  of  this  report. 


FREE  BOUNDARY  PROBLEM  FOR  UNSTEADY  SLAG  FLOW  IN  THE  HEARTH 


Hideo  *>warada* 

1.  INTRODUCTION 

The  hearth  drainage  is  one  of  the  most  important  factors  for  successful 
blast  furnace  operation.  The  slag  is  considered  to  be  more  difficult  to  drain 
than  the  metal  because  of  its  higher  viscosity.  When  the  slag  surface  reaches 
the  level  of  tap  hole,  the  furnace  gas  starts  to  blow  out.  Then  tapping 
should  be  stopped.  The  amount  of  undrained  molten  material  at  the  end  of 
tapping  is  estimated  by  the  shape  of  the  slag  surface.  In  order  to  determine 
the  influence  of  tapping  conditions  due  to  the  shape  of  the  slag  surface,  the 
three-dimensional  problem  of  the  slag  flow  during  tapping  was  solved  by  using 
the  finite  element  method  by  Ichihara  and  Fukutake  [2] .  They  concluded  that 
their  computation  scheme  is  not  efficient  in  practical  use.  This 
computational  instability  was  resolved  by  Kawarada  and  Natori  [11],  using  the 
penalty  method  developed  by  themselves  [5-7,  10]. 

The  objective  of  this  report  is  to  give  mathematical  justification  of 
penalty  formulation,  i.e.,  to  prove  the  convergence  of  the  penalized  free 
boundary  to  the  one  of  an  original  problem  when  we  let  the  penalizing 
parameter  e  tend  to  zero.  In  section  2,  we  review  the  formulation  for  two- 
dimensional  problems  of  the  slag  flow.  In  section  3,  we  give  the  penalized 
formulation  by  using  the  method  of  integrated  penalty.  Section  4  is  devoted 
to  the  assumptions  and  the  main  theorem.  In  section  5,  we  prepare  some 
propositions  needed  to  prove  the  main  theorem,  in  which  the  perturbation  of 
the  penalized  solution  with  respect  to  the  domain  is  discussed.  Finally,  we 
give  the  proof  of  the  main  theorem  in  section  6. 
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2.  FORMULATION 


We  consider  two-dimensional  problem  of  the  slag  flow  in  the 
is  bounded  by  impermeable  boundaries  y  =  0,  x  =  0  and  x  =  a 
1).  One  of  the  vertical  boundaries,  x  =  0,  has  a  tapping  hole 
bottom,  which  we  denote  by 


hearth,  which 
(c.f.  Figure 
near  the 


Figure  1 


We  assume  that  the  hearth  is  packed  with  coke,  through  the  bed  of  which 
the  slag  flows.  Then  Darcy's  law  can  be  applied  for  the  flow  of  slag: 

(2.1)  V  -  -dV* 

where  V  denotes  velocity  of  slag.  The  potential  $  is  defined  as  follows: 


(2.2) 


* 


y 


where  P : 

g: 

P* 


density  of  slag, 
gravitational  acceleration, 
pressure  of  slag, 


2- 


Pgj  pressure  at  reference  point, 
y:  vertical  height  from  a  horizontal  plane, 
d:  permeability  of  slag. 

If  we  substitute  V  into  the  equation  of  continuity: 

(2.3)  div  V  -  0  , 
then  we  have 

(2.4)  A*  -  0  in  the  slag  region  ft. (t)  (t  e  [0,T))  , 

9 

under  the  condition  d  *  constant  and  T  is  tapping  period,  ft, (t)  is 

9 

defined  as  follows: 

(2.5)  ft^(t)  •  {(x,y)|0  <  x  <  a,  0<y<  9(x,t)} 

where  y  »  9(x,t)  represents  the  height  of  the  slag  surface,  which  is  a  free 
boundary,  the  boundary  conditions  for  the  potential  4  are  given: 


(2.6) 

9-0 

y 

on 

y  «  0  , 

(2.7) 

*  «  0 
x 

on 

x  »  0  and  x  *  a,  except  on  the  tapping  hole  VQ  , 

(2.8) 

*  -  V  ^ 

x  out 

on 

ro  ' 

(2.9) 

•  **  y 

on 

y  *  4(x,t)  . 

The  drainage  rate 

Vout 

of  a  working  blast  furnace  increases  as  the  tap  hole 

is  eroded  during  tapping: 

V  _  -  V  (0  <  t  <  T) 

where  «  average  drainage  rate, 

k:  tap  hole  opening  rate, 
t:  tap  hole  erosion  rate. 


(2.11) 


If  -  -VI  ■  -i»  - 

!y-9(x,t) 


+  v. 


-  -d»/  1  +  ♦  +  V, 

x  3n  ,,  ..  in 

y-9(x,t 


-3- 


f 


Where  n  is  outward  normal  to  fl^(t)  and  V^n  ia  the  inflow  velocity  of 
slag  at  the  surface.  The  initial  shape  of  the  free  surface  is  given: 

(2.12)  4»<x,0)  -  *Q(x)  (0  <  x  <  a)  . 

In  this  report,  we  deal  with  the  case  VQUt  is  independent  of  t,  but 
dependent  of  y.  Hereafter,  we  denote  the  free  boundary  problem  (2.4)-(2.12) 
by  (P).  By  using  the  relation  (2.2),  (P)  is  reformulated  into  the 


equations 

for 

{p  and  y  » 

: 

(2.13) 

Ap  * 

0 

in  fl.(t)  (0  < 

t  < 

T)  , 

(2.14) 

fc- 

-P9 

on  y  *  0  , 

(2.15) 

&- 

0 

on  x  «  0  and  x 

*  a 

except  on 

(2.16) 

fe- 

Pg*Vout 

°n  r0  , 

(2.17) 

P  “  i 

po 

on  y  *  $(x,t)  , 

(2.18) 

li¬ 

st 

-d* ( 1  +  / 

pg 

1  +  *1  *  §5|  )  + 

ly** 

Vin 

i 

(2.19) 

<Mx, 

0)  -  *0(x) 

(0  <  x  <  a)  . 

For  simplicity,  we  take  pg  *  1  and  pQ  •  0  and  denote  the  above  problem 


AM  APPROXIMATION  OF  (P»)  BY  MEANS  Of  THE  METHOD  OF  INTEGRATED  PENALTY 


When  we  try  to  solve  (P'),  the  numerical  procedure  must  contain  a 
routine  for  solving  the  potential  problem  (2. 13)- (2. 17)  for  a  given  free 
boundary  y  “  4(x,t).  After  this  is  done,  the  outward  normal  derivative  of 
the  potential  on  the  free  boundary  can  be  calculated.  And  then,  by  solving 
(2.18)  and  (2.19),  the  subsequent  shape  of  the  free  boundary  is  obtained  and 
so  on.  If  we  apply  the  method  of  integrated  penalty  to  solve  the  potential 
problem  (2. 13)- (2. 14) ,  then  the  outward  normal  derivative  of  the  potential 
function  on  the  free  boundary  are  easily  approximated  [3,  10].  This  is  the 
reason  why  we  apply  the  method  of  integrated  penalty  to  free  boundary 
problems . 

Let  B  *  {(x,y)|0  <  x  <  a,  0<y<b}  and  Y(t)  (ten)  be  the 
heaviside  function.  Then  we  penalize  (P* )  as  follows:  Find 
Ip  and  y  “  ♦  }  for  V  e  >  0, 


(3.1) 

A  e  1 

Ap  -  -  . 

Y(y  - 

*€<x,t))pe  -  0 

in 

B  , 

(3.2) 

3pE 

aT-  0 

on 

x  *  0  and  x  = 

a. 

except  on  rQ  , 

(3.3) 

3pe 

3x™  *  Vout  °" 

ro  ' 

(3.4) 

|p!.  .1 

3y 

on 

y  -  0  , 

(3.5) 

pC  *  0 

on 

y  ■  b  . 

Obviously, 

*<y 

-  ♦E(x,t)) 

(0  <  x 

<  a,  0  <  t  <  T) 

is 

the  characteristic 

function  of 

B  - 

ft  -it),  where 

* 


0  c(t)  *  { (x,y) |0  <x<a,  0<y<  $E(x,t)}  . 

♦ 

It  is  well  known  that  if  we  let  e  be  small  enough  in  (3.1),  then  pE 

e 

approximates  p  in  0  (t)  and  p  is  nearly  equal  to  zero  in  B  -  (2  (t) 

♦  ♦* 


V 


[9].  Therefore  the  boundary  condition  (2.17)  is  approximately  satisfied.  If 
we  use  the  method  of  the  integrated  penalty,  the  equation  (2.18)  is 
approximated  by 


(3.6)  *  -d{i  -  r  /  pe(x,n)dn}  +  v 

This  approximation  is  based  on  the  following  fact;  put 


(3.7) 


~  *  V(y  -  ♦E)pE 


(3.8) 


e 


s 


b 

/  re(x,n,t)dn  . 


y 

By  an  application  of  Theorems  1.1  and  1.2  in  [3],  we  have 


(3.9)  r£  -  A  +  *2  •  I  in  V(B)  , 

X  3n|y-t(x,t) 


(3.10)  s£  +  -  /  1  +  •  ( 1  -  Y)  in  P*  (B)  as  e  +  0  . 

x  3n  y=*(x,t) 


By  using  (3.10),  we  have 


(3.ii)  -  /i  +  $1  •  If  I  *  se(x,<)»e)  -  7  /  pe(x,n)dn 

‘,n{y=*(x,t) 

in  a  suitable  topology.  Substituting  (3.11)  into  (2.18),  we  have  (3.6).  Now, 

e  e  e 

b  (x,$  )  is  called  the  integrated  penalty.  Hereafter  we  denote  by  (P  ) 


the  penalized  problem  defined  above.  Finally  we  should  note  that  there  holds 


the  same  discussion  about  (P)  as  mentioned  above 


4.  THE  ASSUMPTIONS  AND  THE  MAIM  RESULT 

4.1.  In  order  to  obtain  our  main  result,  we  need  the  following  assumptions: 
Assumption  (A) 


There 

exists  a  unique  solution 

{p  and  y  ■ 

♦(x,t)} 

of  (P')  such  that 

(4.1) 

p  e  C3'“(fl^(t)) 

for 

t  e  (0,T] 

and  a 

e  (0,1)  , 

(4.2) 

♦  e  C3,°(0,a) 

for 

t  e  (o,T] 

and  a 

e  (0,1)  , 

(4.3) 

♦  e  c^o.t) 

for 

x  e  (0,a) 

e 

Naturally,  vout  included  in  the  boundary  condition  should  have  regularity 
property  in  order  to  obtain  (4.1). 

Assumption  (B) 

There  exists  a  unique  solution  {pe  and  y  -  $e(x,t)}  of  (pe)  such  that 


(4.4) 

pe  e  C3,a(0  ft)) 

for  t  e  (0,T]  and  u  e  (0,1)  , 

(4.5) 

♦e  e  c3,a(0,a) 

for  a  e  (0,1)  , 

(4.6) 

♦e  e  c1(0,T) 

for  x  e  (0,a)  , 

moreover 

(4.7) 

*♦*1  ,  <  X 

C3,°I0,.>  - 

where  K 

is  independent  of  e. 

4.2.  Then  we  have  the  main  result. 

Theorem,  under  the  assumptions  (A)  and  (B),  let  c  ♦  0,  then 

♦  ♦  uniformly  in  (Q,a)  *  (0,T]. 


(4.8) 


5.  PREPARATIONS 


In  order  to  prove  our  main  theorem,  we  have  to  prepare  three 
propositions . 

5.1.  He  consider  the  following  penalized  problem  defined  in  B:  For 


(5.1) 

“APE  +  ~  **<y  ' 

■  4><x))* 

e 

p  «  0  in  B  , 

(5.2) 

3PC 

ST  "  0 

on 

x  -  0 

and  x  ■  a,  except  on 

,  e 

(5.3) 

v 

ox  out 

on 

ro  ' 

.  e 

(5.4) 

-1 

3y 

on 

y  -  0 

9 

(5.5) 

pe  *  o 

on 

y  =  b 

e 

Where  y 

»  $  (x)  is  smooth  in 

i  (0,a) 

and  satisfies 

(5.6) 

0 

<  $(x) 

<  b  (0  <  x  <  a)  . 

Define 

-  {( 

x,y) |0 

<x<a,  0<y<  $ ( x) } 

and 

r 

-  { (x,y) | 0 

<  x  <  a,  y  -  ♦(x)}  . 

Let  p° 

be  the  solution 

of  the  problem: 

(5.7) 

Ap°  ■  0 

in 

• 

(5.8) 

a  0 

fr-* 

on 

x  -  0 

and  x  *  a,  except  on 

(5.9) 

£.v 

3x  out 

on 

ro  ' 

(5.10) 

p°  -  0 

on 

r  and 

3° 

■  -1  on  y  -  o  • 

It  is  well  known  [9]  that 


-8- 


(5.11) 


pE  ♦  p^» ( 1  -  Y)  strongly  in  H^(B)  as  €  +  0 


Then  we  have 

Proposition  1 .  Let  e  be  small  enough  and  m  >  0.  Then 


<5-’2>  lpE,n,r  *  <’</e> 

<5-131  ,p‘  "  /r  '  f5-'.,r  i  0le’  ‘If-'-.r  ■ 


(5.14) 


i/  Pedri 


Pe,m,r  ^  0<e),Pe,m.i,r  • 


Proof.  See  Kawarada  and  Hanada  ([4],  Chapter  2,  p.  4). 

5.2.  We  fix  some  y  “  ♦°(x),  which  satisfies  (5.6)  and  is  smooth  enough  and 
we  denote  the  problem  (5.7)-(5.10)  by  Pp.  Then  any  domain  ft^  which  is 
diffeomorphic  to  ft  Q  can  be  defined  enough  near  ft  Q  such  that 

♦  t 

(5.15)  0^  -  t(x,y)|0  <x<a,  0<y<  *(x)} 

where  $(x)  *  $°(x)  +  6$(x)  (0  <  x  <  a) 

with  6$  e  CW(0,a)  (p  >  1)  . 


We  consider  a  family  of  the  penalized  problem  (5.1)- (5.5),  in  which 

Y(y  -  $(x))  is  replaced  by  Y(y  -  $'(x))  ($’■$+  6$).  Denote  the 

e* 

solution  of  this  problem  by  P  .  Then 
Proposition  2. 


(5.16) 


lim 

1641  >0 

CM(0,a) 


exists  , 


£  ft 

(5.17)  e  C  (B)  (a  e  (0,1))  . 

(*)  I VI  _  indicates  the  norm  of  V  in  H^G) . 
fll^G 


r 


Proof .  See  Dervieux  (II]/  Chapter  3,  p.  47) , 
£ 

properties  as  p  does  in  Proposition  t . 
Proposition  3. 


Also 


3$ 


satisfies  the  same 


(5.18)  «|J-  ~  ~  Pe<n,ey>lm  r  <  0(/T)l<n.ey>lm+1  x 

(5-19)  ■/,  it  dpi  -  ^  ff-'m.r  ^  0(e),<n'V'*+i,-,r  x 

r*  A 


where  n  is  unit  outward  normal  to  £).,e  *  Vy  *  ( 0, 1 ) ,  <•  ,  •>  implies  the 

#  y 

indicates  the  norm  of  v  in  <n. 


inner  product  in  and  lvlf 

e 

Proof,  p  satisfies 
(5.20) 

Differentiating  both  sides  of  (5.20)  with  respect  to  4>,  we  have 


m,-,r 


-ApE  +  ~  *  Y(y  “  ♦(x))p£  *  0  in  B  . 


(5.21) 


'A  H"  +  e  Y(y  "  ♦>  dt  ~  e  *  I?  ’  p£  =  °  in  ® 


F  3n  1 

Put  Q  =  t*?-,  multiply  ♦  6  H_(B)  both  sides  of  (5.21)  and  integrate  in 
3<p  0 


B.  Then  we  have 


(5.22) 


/{  VQe*V*  dxdy  +  7  //  YQ%  dxdy  -  7  //  4“  P%  dxdy  «  0 


*  "  8y 


By  using  the  relation; 


(5.23) 


//  f“  P£*  dxdy  -  /  pb<n,ey>*  dr  , 


we  have 


(5.24)  Implies 


(by  Green's  formula)  where 


and 


(5.25) 

(5.26) 

(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 


AQ0  -  0  in  , 


-AQ®  ♦  7  Q®  -  0  in  B  -  Q*  , 


Oq  “  °1  on  F 


3 oc  3oe 

w0  U1  1  e  v 

-  “  "  7  *  P  o»  f 


3n  3n  e 


!S 

3x 


on  x  »  0  and  x  -  a,  except  on  1*.  (i  *  0,1)  , 


3V 


out 


3x  3y 

3^ 


on  rfl  , 


“0  on  y  -  0  , 


Q®  «  0 


on  y  •  b 


(5.27)  comes  from  (5.17).  If  we  use  the  method  developed  in  Kawarada  and 
Hanada  ((4],  Chapter  5,  p.  15),  we  are  able  to  obtain  (5.18)  and  (5.19). 


6.  THE  PROOF  OF  THE  MAIN  THEOREM 

6.1.  He  reformulate  (2.18)  into  the  following  penalty  formulation  by  means  of 
the  method  of  integrated  penalty: 


(6.1) 


li  s  -<j (i  ~  j  J  qe(x,n)dn)  +  v4w  +  0{/e) 


in 


q  is  the  solution  of  the  problem: 

(6.2) 


-AqE  +  j  •  Y(y  -  $(x,t))»qC  “  0  in  B 


£  £ 
and  q  satisfies  the  same  boundary  conditions  as  p  does  on  9B.  In  order 

to  justify  (6.1),  we  may  show  the  following  inequality: 


qe<*  +  S'  ♦  ♦*  *  f?,C(D 


<  l/  1  +  *2I  *1“  /  q^dT1  +  |EI< 

*  x  c(0,a)  e  .  3n  < 

r 

<  l/  1  +  *2I  .  •  I—  /  qedTi  +  |£|  . 

-  Tx  C(0,a)  e  3n  H^( 

<  ./  1  •<*>  (by  (5.U)  and  (5..4H 


C(D 


'(D 


H  (D 


<  o(/?>  ♦  ♦iic,0,.)-u>cl(r( 

<  oc/o  ./>  M|ic(0ill'ip> 


(t  e  (0,t) )  . 


C  (fl. (t)) 


If  we  use  the  regularity  property  of  the  solution  of  an  elliptic  boundary 
value  problem  (see  Theorem  3.3.1  in  [8]),  we  have 


(6.4) 


'pl  3  0  <  C(l*l 

C  '  (fl^(t))  C3'°(0,a) 


)  •  HV  I  _ 

out  c2.«(ro, 
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Combining  (6.3)  and  (6.4)/  we  have 


(6.5)  |i/  ,*«  •  |£|c(r)  <  0(/T).C(l6lc3,a(0  >))  <  0(/<) 

♦ 


Here  we  used  (4.7)  in  the  assumption  (B) . 
6.2.  $E  satisfies 


(3.6) 


fi_  .  1  _  1  J  pe(x,n)dn)  +  vin  . 


e  e 

Substracting  (6.1)  from  (3.6)  and  putting  M  “  i  -  $ »  we  have 


(6.6) 


f“"  -  -  ~  <F(*6)  -  ?(♦))  ♦  0(/e) 


b  b 

where  F(4€)  -  /  pCdn  and  F($)  »  /  qEdh  .  Since  there  holds 


F(*E)  -  F(*)  -  /  U  (s*E  +  (1  -  s)*)ds‘ME  , 
0  d* 


we  may  estimate  the  following  FrSchet  derivative: 


b  „  e 


(6.7) 


dF  ,  |p_  ef ,  e. 

«  ml.  w  *  ' p  ’  ' 


e  e 

where  ♦  *  a$  ♦  (1  -  s)$  (s  e  (0,1)).  From  the  assumptions  (A)  and  (B),  we 


(6.8) 


«1>  I  ,  „  <  K  (a  e  (0,1))  . 

C3,a(0,a) 


bet  us  recall  that  tt-  was  defined  in  proposition  3  of  section  5. 
c  3p  e 

6.3.  If  we  put  Q  *  then  Q  satisfies  some  properties  included  in 
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propositions  2  and  3.  Define  T  »  r  (t)  *  i(x,y)|0  <  x  <  a,  y  ■  +  (x,t)}. 
Let  us  estimate  (6.7)  for  m  >  0, 


(6.9)  if  QCdT)  -  pe(<J>e)l 


Hm(re) 


<  */  i  +  •(/  Qea(re)i  -  /e  Q£) 

(rV 


♦  /e  /  1  +  •  (Q6  -  —  <n,e  >p€)l 

x  /e  y  HB(re) 


(note  /  1  ♦  ^  •  <n,e  >  «  1) 

x  y 


<  "  1  +  ♦*  1  «  -  e  •('/  0Ed(rV  -  /e  c£» 
w  <r  >  (rV 


in  e 
Hm(r) 


+  /e  •  IQE - -  pE<n,e  >1  _  ) 

/c  y  h  (r  ) 


<  0(e)  l/  1  +  4»e  I  •  l<n,e  >1  ..  *I~E| 

x  wm,"(r  )  y  W**1'*!'6) 


(by  (5.18)  and  (5.19)).  Putting  m  «  1  in  the  above  inequality. 


(6.10)  if  Qedn  -  Pe«pE)i  <  e-C(l*€l  )*l|fl 

C(T  )  w3'*(0,a)  3  H2 (r  ) 


<  e»c(l*El  )  <  eK  . 

<T'  <0,a) 


Here  we  used  the  same  discussions  as  in  6.1  and  note  that  K  is  independent 


of  s  e  (0, 1 ) 


6.4.  Integrating  (6.6)  in  (0,t),  we  have 


(6.11) 


M6  -  -  7  /  (r(#e)  -  r(#))dr  ♦  0(/7)»t 

£  0 


tram  which 


a  r*  .a 


6.12)  ■«*<•. if/ 


Substituting  (6.10)  into  (6.12),  we  have 


(6.13)  '"‘<••‘"0(0,.)  S  d,‘  l  *  0,/'1 


from  which 

(6.14)  Hl!e(«,t),c(ofa)  i  0(/*>  for  vt  e  l0»T] 


Therefore 


(6.15) 


teio,T] 


'*  '-^'ao,.)  *  0  a*  E  *  0  * 
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